Optical parametric oscillators in isotropic photonic crystals and cavities: 

3D time domain analysis 
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We investigate optical parametric oscillations through four-wave mixing in resonant cavities and 
photonic crystals. The theoretical analysis underlines the relevant features of the phenomenon and 
the role of the density of states. Using fully vectorial 3D time-domain simulations, including both 
dispersion and nonlinear polarization, for the first time we address this process in a face centered 
cubic lattice and in a photonic crystal slab. The results lead the way to the development of novel 
parametric sources in isotropic media. 



O 
• I— I 

> 

Oh 

q , 

o ■ 

^ ■ 



> 

o : 
m . 
o ■ 

^ : 
o , 

o . 

^ ; 

Oh 



Photonic crystals (PC) are synthetic structures de- 
signed for specific optical applications, realized by a (2 or 
3D) periodic pattern embedded in a bulk material. PC 
exhibit forbidden and allowed frequency bands of states 
for electromagnetic waves, depending on dielectric con- 
stants and geometric features, light polarization and di- 
rection of propagation. The analogy to solid state crys- 
tals in electronics is at the origin of their name 
An alternative way to view PC is in terms of a peri- 
odic tight distribution of resonant cavities (RC), such 
that the resulting system presents bands of states owing 
to mode splitting in the elemental RC. A description of 
these structures exclusively in terms of free propagation 
within the allowed bands is rather limiting, especially in 
high-index contrast and finite-dimension PC. 

The resonant response of PC is manifestly relevant 
when considering nonlinear optical processes. By their 
very nature, in fact, the latter can couple cavity modes, 
yielding a spectrum which is expected to significantly 
depend on the specific distribution of states. If sufficient 
energy is available, the interplay between frequency (or 
mode) mixing and cavity effects may build up oscillations 
driven by the nonlinear gain. In PC, a high Q-factor |^ 
entails low thresholds and, hence, a broad output spec- 
trum which peaks at the most state-crowded frequencies. 

While quadratic nonlinearities are the basis of most 
optical parametric oscillators (OPO), the process of four- 
wave mixing (FWM) in a cubic medium can trigger non- 
linear gain and hence oscillations, as well, (see e.g. 
^, 7, 8, 9], and references therein) Since it is present in all 
dielectrics, a cubic nonlinearity is an excellent candidate 
for highly integrated parametric sources in PC, their re- 
alization crucially depending on materials. Furthermore, 
the large density-of-states (DOS) characteristic of PC re- 
laxes the selection rules that govern FWM parametric 
gain in resonant structures, |a| while providing several 
synchronized (phase-matched) interactions. The strict 
conditions that affect FWM in microresonator (without 
PC) have been pointed out in 10], where the observa- 
tion of very low threshold (due to high Q-factor) Raman 
lasing has been reported. 



In this Letter we present a general model of OPO 
through FWM in photonic crystals and cavities real- 
ized in isotropic materials (and hence lacking a quadratic 
nonlinearity). Such possibility is then verified by means 
of extensive vectorial time domain simulations, using a 
finitc-diffcrence-time-domain (FDTD) parallel algorithm 
which accounts for material dispersion in both the linear 
and nonlinear regimes. To the best of our knowledge, 
these are the first fully comprehensive 3D vectorial time- 
domain simulations of cavities and PC encompassing dis- 
persion and a cubic optical response. 

Nonlinear optical phenomena in high index contrast 
3D-PC were previously investigated numerically with ref- 
erence to an instantaneous third-order nonlinearity and 
to its role in distorting the band-structure 0, 0| . A 
general theory was developed in using the "effec- 
tive medium" approach but leaving aside the resonant 
features. More recentl y, s ome experimental studies have 
been reported in PC [3| (see also 0| for an updated 
review on nonlinear photonic crystals). 

Hereby we focus, both theoretically and numerically, 
on frequency coupling associated to collective phenom- 
ena between several states in cavities. Let us consider 
an RC, or a PC of finite extent, realized in an instanta- 
neous isotropic medium. The PC is defined by a position- 
dependent refractive index n(r), in a region of volume V 
and surface S over which periodic boundary conditions 
are enforced. Electric and magnetic fields can therefore 
be expandend in terms of cavity eigenfunctions, denoted 
by a generalized index /i, according to 

E(r, t) = J2f, af,{t)'Ef,{r) exp{~iuj^t) 
H(r, t) = a^(OH^(i") exp(-itj^i). 

In as usual, the summation is extended over nega- 
tive frequencies in order to make the fields real-valued, 
with ~fi corresponding to complex conjugation. Apply- 
ing the Lorentz reciprocity theorem to an instantaneously 
respondin g m edium, one can derive the coupled mode 
equations [l^ for the amplitudes a: 



(1) 



dt 



■exp(iw^t)(E* • J) 



(2) 
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where the brackets denote the integral over V", and J is 
the current density associated to the nonhnear polariza- 
tion Pnl (i- e. J = 9tPjvi) or other coupling mecha- 
nisms, if present. In Q we also included a loss coeffi- 
cient Ku., such that the quality factor is — ct;^/2K^. 
illllll 

For a third-order material we have: Pnl ^ 
x'™"°(r)i?™iJ"i?°, the superscripts denoting the carte- 
sian components, and having omitted the sum over re- 
peated symbols. Eqs. (O then become (the dot stands 
for time derivative): 



dt 

a/37 



(3) 

with coupling coefficients defined as the overlap integral 
between Bloch functions and the nonlinear response, i.e. 

Being interested in optical parametric oscillations due 
to a cubic nonlinearity, we consider a mode of ampli- 
tude a and frequency w which, pumped by some external 
source, couples with two other modes of amplitues a_|_ and 
a_ at frequencies uj± respectively, such that cl;_|_ + cj_ = 
2u!. Therefore, in the stationary regime © reads (for 
simplicity we neglect self- and cross-modulation terms): 



P 



(4) 



with g the relevant nonlinear coefficients and p the pump- 
ing mechanism (see 0|). The Mauley- Rowe relations are 
P+/uj+ = P-/uj- = -P/2uj + pa* , with P = 2K|ap the 
output power corresponding to field amplitude a (and 
similarly for P± ) . "Below threshold" solutions of Q are 
a± = and a = p/n (fixing arg(a)), while the threshold 
for oscillation is 
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uj+u- \g\'' 



(5) 



or, in terms of the output photon flux at w: P/lj — 
l/2\g\Q^jQ+Q^. As expected, large Q's reduce the 
threshold, while abobe threshold the value of a is 
clamped to (O and 



l«±l 




(6) 



As in standard OPO's, the latter shows that the excess 
pumping (Ap = pj na— 1) is transferred to u;±. 

Once the system is oscillating all the modes will vibrate 
at frequencies originated by four-wave mixing of uj and 
uj±. Henceforth, in the no-depletion approximation, a 
and a± can be treated as source terms in 0. If ijJ± — 
uj ± Aw, the generated frequencies are Sw, ± Aw, 
uj ± 2 Aw, 2lo ± Aw. Thereby, a^, is the sum over all the 



FWM terms at each w^. For instance, if uji = w+w+— w^, 
for any mode of order v we have 

+ Ai^fly = iW2/t,(w, +W+, -w_)aala+e*^"''~'^^''*, (7) 

being fu{uj, +^+- — w_) an effective nonlinear coefficient. 
The stationary solution of Q is 



Ki, + i{uJu - W2) 



(8) 



The energy at W2 is £{uJ2) = Yl,u \onf^\'^i which can be 
cast in the form 



w||/.(w2)|LgAp^ 
32|5|3Q2yQ7Q: 



p(w) 



K2 -t- (w — W2)' 



;dLO, (9) 



where we introduced an averaged nonlinearity |/i/(w2)|avg 
and the DOS p(w) of the cavity. In deriving ^ we 
considered that, if Q{lo2) = Q2 is sufficiently high, the 
relevant modes are those in the proximity of W2. If the 
DOS is smoothly varying (see [13), expression states 
that the energy is approximately proportional to p(w2). 

Increasing the pump fluence, the amplitudes a± be- 
come sources of other parametric processes, and more 
frequencies are generated. At high powers the whole os- 
cillation spectrum will resemble the DOS (with the obvi- 
ous exception of those interactions for which g = due 
to symmetry). It will be peaked in regions where the 
states are denser, i.e. in the proximity of the photonic 
band-gap (PBG) for a PC or near cut-off for an RC.[l^ 

It is important to determine the bandwidth of the 
pump modes. The pump, at frequency w^, transfers 
energy to each mode at w = Wi„ — (5, according to its 
Lorcntzian lineshape. We can account for this effect by 
writing the source term as pe~'^^* and solving the time- 
dependent coupled equations. This leads to replacing k 
and LJ hy K + iS and w^, respectively, in the first of 
and in corresponding ansatz's for a±. As a result, when 
S ^ Q the threshold value for a increases, and the pump- 
bandwidth is determined by 



\p/^\ 



1 + (S/n) 



(52 



[(«+ + «:_)/2]2' 



(10) 



with Lu± = w+2(5k±/(k_(-+k-). For small 6, with k± ^ k, 
from H10|l we have \S\/uj = ^/Ap/Q. The total pump 
energy is the sum of the energies in each pump-mode: 
the oscillator will be more efficient if the DOS is strongly 
peaked around ujm ■ 

With the aim of demonstrating the feasibility of OPO's 
in isotropic PC/RC, we performed fully vectorial numer- 
ical simulations of nonlinear Maxwell equations, with a 
code based on the FDTD approach for dispersive mate- 
rials |2fij |. Maxwell equations in vacuum were coupled 
with a nonlinear Lorentz oscillator (see Eqs. Hll|) be- 
low), yielding the induced polarization P {P^ = P • P) 
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in regions where material is present: 

V X E = -fiodtU 

V X U = eodtB + dtP (11) 
d^P + 2jodtP + iolf{P)P = eo(e, - l)c^gE. 

In the algorithm, the Yee's grid 0| was used to enforce 
continuity between different media, and uniaxial phase- 
matched layers (UPML) were adopted at the bound- 
aries ^23] (details will be provided elsewhere). Partic- 
ular care is required in adopting a specific form of the 
Lorentz oscillator f{P), with /(P) = 1 describing a 
linear single-pole dispersive medium. For an isotropic 
Kerr material we followed the suggestion by Koga [23| . 
i.e. /(P) = [1 + (P/Po)2]"3/2 is a measure of 
the nonlinearity) . For small P/Pq the latter reproduces 
a Kerr response while also accounting for higher order 
terms. Compared to the standard Kerr /(P) = 1-1- xP^, 
the resulting algorithm is stable even near the Courant 
limit. For the effective Kerr-law coefficient n2 we 
chose 712 = 1-5 X 10~^'^rn^ /W, representative of an en- 
tire class of semiconductors |2J| . This value gives Po = 1 
(MKS units). Because of the significant computational 
resources needed for such a numerical approach, we par- 
allelized the algorithm. (3QJ 

A proper test-bed for our theoretical analysis is pro- 
vided by a dispersive medium with non-instantaneous lin- 
ear and nonlinear responses, modeled by pi(l and solved 
numerically with no- approximations. We computed the 
response of a Face-Centered-Cubic (FCC) lattice (of pe- 
riod A) of air-spheres (radius r = 0.3535A) embedded 
in a dielectric. This is one of the simplest structures 
admitting a complete PBG jl^. The PC, of dimen- 
sions 8/im X 8/iTO X 8/iTO, was placed in air and excited 
by a 2/iTO-waist linearly y-polarized gaussian beam, ob- 
tained by a total field/scattered field laver.|2^ The input 
temporal profile, for the quasi-cw excitation mentioned 
below, exhibited smooth transitions from zero to con- 
tinuous wave, as mimicked by an mnm pulse (m = 2, 
n = lO''). |2y| Its spectrum was well peaked about the 
carrier frequency. The FCC lattice for a material with 
index 3.5, such as Si or GaAs, has a complete band gap 
around a normalized frequency A/A = 0.8, with A the 
wavelength. To get a gap near A = ISOOrtm we chose 
A = 1200nTO, and the parameters of the single pole dis- 
persion were taken as — 11.971392, luq ~ 1.1 x 10^^ 
and 7o = 2 X 10^ (in MKS units), yielding an index = 3.5 
at A = 1500nm. The integration domain was discretized 
with dx dy ^ dz ^ 30nm and time step dt = 0.03/s, 
allowing more than 40 points at each wavelength around 
A = 1500nm and runs with 20000 steps in time (the spec- 
tral resolution is of the order of lOnm at A = 1500nm). 

The DOS is typically calculated by the use of plane- 
wave expansion, omitting material dispersion and for in- 
finitely extended structures (see, e.g., jl^ 1131 )■ Hence, 
it is all but straightforward the application of the stan- 
dard approach to the case under consideration. For this 




wavelength X (nm) 

FIG. 1: (Colors) Spectrum of the transmitted signal for a 
single-cycle low-power excitation. The dotted line is the input 
spectrum (normalized to unity). The star marks the wave- 
length used to pump the OPO. A sketch of the FCC lattice 
is on the upper right. 



reason, we resorted to a time-domain approach to deter- 
mine the states of the FCC-PC. A very low-power (InVF) 
single-cycle pulse excited the PC (along the TX di- 
rection) and the transmitted signal (its i?^ component) 
was analyzed just after it. The resulting spectrum is 
shown in figure ^ where the peaks correspond to con- 
centrations of states (taken aside the low-frequency os- 
cillations due to the finitess of the structure (J|) compat- 
ibly with symmetry constraints. The band structure of 
this medium encompasses a PBG around 1500nm and a 
pseudo-gap around 2400nm. 

In figure |21 we display the spectrum obtained in 
a low-symmetry point at the center of the PC, for a 
beam propagating along the direction TX. Different in- 
put powers were injected with a 600/s quasi-cw excita- 
tion. The pump wavelength A = 1336nrn corresponds to 
A/A = 0.898, i.e. close to a state by the (frequency) up- 
per edge of the PBG, as indicated by the star in figure^ 
As visible in the insets, large output spectra are attained. 
No oscillations appear at frequencies within the PBG, 
with a smooth spectrum in the large wavelength region 
and several peaks above the PBG upper edge. Each peak 
corresponds to a region dense of states, and hence to 
an efficiently generated frequency. 

As a second representative case, we considered a pho- 
tonic crystal slab with a triangular lattice (A = 450nm, 
hole radius r = 135nm), as experimentally investigated 
in l2i|. The slab is 270nm thick, with in-plane size 
2fim X 2/im. It is pumped by a 0.5/iTO-waist TE-like 
polarized gaussian beam incident along the TM direc- 
tion. The pump wavelengths are in proximity of the 
upper and the lower edges of the guided-modes gap (lo- 
cated in the interval 1260-1750nm for the infinite non- 
dispersive structure), and the material constants are 
LOO = 1.1 X 10^^, 70 = 0, Po = 1, = 10.9379, 
and refractive index n = 3.3 as in AlQ,iGao,gAs at 
A — 1800nm. jl^ The discretization was implemented 
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FIG. 2: Oscillation spectrum inside the FCC for Pin = IMW. 
Several excited wavelengths correspond to the spectral peaks 
in figure Q The asterisk marks the pump wavelength. The 
insets are close-ups of the spectral interval around 1200nm 
for different input levels. At high powers a slight shift of the 
resonance peaks is due to self- and cross-phase modulation. 
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FIG. 3: (Colors) Oscillation spectra inside a photonic crys- 
tal slab (upper right), with generation of several wavelengths 
{Pin ~ lOOfeM^). The thick line corresponds to a pump 
at A = 1250nm, by the PBG upper edge; the thin line to 
A = 2000nm, by the lower edge. The inset shows the output 
spectrum for the low power (InW) single cycle excitation, the 
stars indicate the OPO pumps. 



with dx = 15nm, dy = 25nm, dz = 19nm, dt = 0.02/s, 
and both single cycle and cw feeding were carried out for 
up to 30000 steps in time (= 600/s). Figure |3| shows 
OPO spectra of the transmitted field for pump wave- 
lengths A = 1250 and 2000nm, respectively (mnm pulse, 
m = 5, n = 600), and input power Pin = lOOkW. Even 
at such low power, the broadening is enhanced owing to 
the spatial confinement afforded by a planar geometry. 

In conclusion, we have analyzed optical parametric os- 
cillations in photonic crystal cavities. Through fully vec- 
torial 3D numerical simulations we have demonstrated 
that regions of significant mode-concentration and high 
Q-factors favor the onset of parametric oscillations due 
to four-wave mixing gain. Moreover, a properly tai- 
lored DOS can favor the generation of specific frequen- 
cies, making nonlinear amplification possible even with 
isotropic (i.e. centrosymmetric) materials in microcav- 
ities. These results pave the way to novel genera- 



tions of specifically tailored parametric sources employ- 
ing isotropic materials in highly integrated geometries. 
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